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Introduction 


The  launching  of  Landsat  I  in  1972  provided  a  new  tool  for  agri- 
cultural inventories  and  analyses.    Multispectral  scanner  data  became 
available  from  Landsat  I  for  extensive  areas  on  a  repetitive  basis 
every  18  days  and  with  the  launching  of  Landsat  II,  every  9  days.  The 
prospect  of  timely  and  relatively  inexpensive  data  encouraged  agricul- 
turists to  experiment  with  Landsat  in  their  inventory  efforts.  This 
study  is  an  example  of  one  such  effort.    The  objective  was  to  inventory 
new  cropland  development  between  1972  and  1976,  a  period  of  high  agri- 
cultural prices,  to  obtain  a  measure  of  cropland  response  to  favorable 
cost-price  relationships.    In  this  report,  the  results  of  the  study 
are  presented  and  analyzed.    The  limitations  of  Landsat  imagery  for 
projects  of  this  nature  are  discussed  and  guidelines  for  future  Land- 
sat land  use  inventory  projects  are  proposed. 

Review  of  Landsat  Projects 
Related  to  Agriculture 

Since  1972,  numerous  projects  have  been  initiated  using  satellite 
imagery  to  monitor  the  agricultural  sector.    The  basic  element  of  Land- 
sat data  is  the  set  of  measurements  or  (signature)  by  the  satellite's 
multispectral  scanner  (MSS)  of  a  .4  hectare  area  of  the  earth's  surface. 

'v 

The  MSS  measures  the -amount  of  radiant  energy  reflected  and/or  emitted 
from  the  earth's  surface  in  various  bands  of  the  electromagnetic  spec- 
trum.   The  Landsat  II  and  Landsat  III  MSS's  have  four  and  five  bands, 
respectively.    Each  band  can  be  analyzed  separately  or  in  combination 
with  other  bands  to  classify  land  uses.    Two  forms  of  data  are  available 
for  analysis — MSS  imagery  or  MSS  digital  data  on  tapes.    Manual  techniques 
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are  used  to  interpret  the  former  and  computer-assisted  techniques  are 
used  to  analyzed  the  digital  data.    This  section  is  not  meant  to  be  a 
comprehensive  review  of  agricultural  Landsat  projects,  but  rather  a 
review  of  the  approach  of  some  representative  studies. 

The  LACIE  project,  the  Large  Area  Crop  Inventory  Experiment, 
was  a  cooperative  project  of  the  U.S.  Department  of  Agriculture 
(USDA),  the  National  Aeronautics  and  Space  Administration  (NASA),  and 
the  National  Oceanic  and  Atmospheric  Administration  (NOAA)  of  the  U.S. 
Department  of  Commerce.    It  was  initiated  to  develop,  test,  and  evaluate 
a  system  for  predicting  area,  yield,  and  production  estimates  for  ma- 
jor wheat  producing  areas  of  the  world  through  the  combined  use  of 
satellite,  weather,  climate,  and  historical  agricultural  data  (LACIE, 
1977).    The  objective  of  LACIE  was  to  estimate  the  production  of  wheat 
on  a  country-by-country  basis.  Aerial  extent  of  wheat  was  derived  from 
a  computer  assisted  analysis  of  Landsat  data  while  yield  was  estimated 
from  statistical  models  which  relate  crop  yield  to  local  meteorological 
conditions,  notably  precipitation  and  temperature. 

LACIE  was  a  three  phase  project  beginning  in  November  1974  and 
ending  in  mid-1978.    During  Phase  1,  the  LACIE  experimental  system  was 
put  together  and  tested  on  the  1974-75  U.S.  Great  Plains  wheat  crop. 
Phase  2  tested  LACIE' s  ability  to  estimate  acreage,  yield,  and  produc- 
tion of  the  1975-76  wheat  crop  in  the  U.S.  Great  Plains,  Canada,  and  two 
indicator  regions  in  the  Soviet  Union.    During  Phase  3,  LACIE  produced 
in-season  acreage,  yield,  and  production  estimates  for  the  U.S.  Great 
Plains  and  the  total  USSR  wheat  crop. 

Evaluations  of  LACIE  through  USDA  (Statistical  Reporting  Service) 
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ground  survey  estimates,  indicate  the  project  was  more  successful  in 
the  winter  wheat  regions  of  the  world  as  there  was  difficulty  in  dif- 
ferentiating between  spring  wheat  and  other  spring-sown  small  grains. 
In  addition,  the  regression  yield  models  formulated  were  not  adequately 
responsive  to  episodic  events  and  require  improvement  to  achieve  accu- 
rate estimates  in  years  with  large  deviations  from  normal. 

The  Statistics  group  of  USDA  has  utilized  Landsat  data  as  an  aux- 
iliary variable  for  its  agricultural  surveys  (Sigman,  197  8).  Crop 
acreages  are  estimated  by  a  regression  estimator  with  enumerator  data 
from  the  June  Enumerative  Survey  as  the  primary  survey  variable  and 
Landsat  digital  data  as  the  auxiliary  variable.    Major  demonstration 
efforts  have  been  conducted  in  Illinois,  Kansas,  and  Kings  County,  Cali- 
fornia using  this  technique.    All  of  these  studies  have  been  research 
oriented,  and  except  for  the  1977  study  in  Kings  County,  California, 
none  of  these  demonstration  projects  have  produced  timely  crop  acreage 
estimates.    Based  on  this  methodology,  however,  Landsat  crop  acreage 
estimates  are  expected  to  be  used  for  final  season  estimates  in  Iowa 
in  1978. 

Landsat  digital  data  has  been  utilized  in  the  state  of  North  Dakota 
in  a  Regional  Environmental  Assessment  Program  (REAP)  to  produce  a  land 
cover  inventory  for  the  state  (REAP,  1977).    Multi-color  land  cover  maps 
for  ten  land  cover  categories  were  the  final  product.    This  project 
found  several  categories  particularly  difficult  to  map.    Fallow  fields 
and  built-up  areas,  for  example,  reflected  identically.    Selected  up- 
dating  of  the  maps  is  planned. 
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The  use  of  Landsat  has  been  especially  investigated  in  less  devel- 
oped countries  where  secondary  ground  data  is  usually  not  available. 
One  example  is  the  use  of  Landsat  data  in  estimating  agricultural  pro- 
duction in  LAOS  (Hooley,  1977).    Manual  interpretation  techniques  were 
used  to  provide  data  that  was  relatively  cheap  and  did  not  require  the 
availability  of  highly  specialized  computer  processing  capabilities. 

These,  and  other  similar  projects,  have  used  Landsat  imagery  to 
estimate  crop  acreage  in  an  inventory  effort  or  in  combination  with 
other  information  to  estimate  crop  production.    Landsat  was  used  as  a 
primary  data  source  in  most  of  the  studies,  although  the  USE&  Statis- 
tics projects  have  used  Landsat  data  as  auxiliary  variables.  Inte- 
grating Landsat  with  other  measurement  techniques  may  be  the  most 
effective  use  of  satellite  imagery,  however  many  projects  seek  infor- 
mation not  available  through  other  methods. 

The  study  presented  in  this  report  uses  Landsat  as  primary  data 
as  did  most  of  the  studies  reviewed.    The  time  period  studied  and  the 
data  sought  did  not  correspond  to  existing  periodic  inventories.  Hew- 
ever,  this  study  differs  in  purpose  and  methodology  from  previous 
studies.    The  purpose  is  to  determine  land  use  change  rather  than  in- 
ventory current  land  use.    Two  time  periods  rather  than  one  are  em- 
ployed.   Landsat  imagery  for  1972  and  1976  was  utilized  to  estimate 
cropland  development  during  four  years  of  favorable  price-cost  rela- 
tionships for  farm  products. 

Description  of  Study 

Study  Design 

Two  regions  were  chosen  for  analysis:  the  Central  Great  Plains 
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Winter  Wheat  and  Range  Region  (Great  Plains  region)  and  the  Southeast- 
ern Alluvium  Coast  and  Tidewater  Region  (Southeastern  region).  These 
regions  were  selected  for  investigation  on  the  basis  of  evidence  of 
cropland  development — irrigation  development  in  the  Great  Plains  and 
clearing  of  forestland  for  crops  in  the  Southeast.    A  systematic 
sample  of  20  percent  of  the  counties  in  each  region  was  chosen  for 
analysis  (table  1).    Forty-five  counties  were  chosen  for  analysis  in 
the  Great  Plains  and  85  counties  were  selected  in  the  Southeast. 

Cloud  cover,  a  problem  inherent  in  Landsat  imagery,  introduced 
some  bias  into  the  selection  process.    Cloud  cover  prevented  the 
inclusion  of  some  first  choice  counties  in  the  sample.  Alternate 
counties  were  chosen  for  13  .percent  of  the  counties  in  the  Great 
Plains  and  19  percent  of  the  counties  in  the  Southeast. 

Land  use  change  between  two  time  periods,  1972  and  1976,  was  in- 
vestigated.   Five  categories  of  land  use  change  in  the  two  regions 
were  to  be  identified.    In  the  Southeast  the  categories  were  conver- 
sions of  (1)  forest  to  cropland  and  (2)  pasture  to  cropland.    In  the 
Great  Plains  the  categories  were  conversions  of  (1)  nonirrigated  crop- 
land to  irrigated  cropland,  (2)  rangeland  to  irrigated  cropland,  and 
(3)  rangeland  to  nonirrigated  cropland. 
Procedures 

The  interpretation  of  the  Landsat  imagery  was  contracted  to  the 
Resource  Information  Laboratory  (RIL)  at  Cornell  University.  Manual 
techniques  developed  by  the  RIL  in  earlier  research  projects  were 
adapted  to  this  project  (Hardy,  1975). 

Landsat  imagery  was  obtained  for  the  study  counties  from  the  EROS 
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Table  1.    State  and  regional  groupings  of  130  sample  counties 


Alabama 


Southeastern  Alluvium  Coast  and 
Tidewater  Region  (85  counties) 

Arkansas 


Barbour  (alternate  for  Bullock  AL) 

Choctaw  (alt.  for  Coffee  AL) 

Dale  (alt.  for  Elmore  AL) 

Henry 

Mobile 

Tuscaloosa 

Wilcox  (alt.  for  Clark  AR) 


Desha 

Cleveland  (alt.  for  Hempstead  AR) 
Lee 

Little  River  (alt.  for  Nevada  AR) 

Phillips 

Prairie 


Florida 

Baker 

Calhoun 

Dixie 

Jefferson 

Liberty 

Nassau 

Washington 


Georgia 

Brantley 

Brooks 

Clay 

Clinch 

Dodge 

Evans 

Houston  (alt. 

Jeff  Davis 

Macon 

Pierce 

Pulaski 

Seminole 

Terrell 

Turner 


for  Jenkins  GA) 


Kentucky 
Calloway 


Louisiana 

Beauregard  (alt.  for  Bossier  LA) 

Calcasieu  (alt.  for  Vermilion  LA) 

Caldwell  (alt.  for  Evangeline  LA) 

Cameron  (alt.  for  Lafourche  LA) 

Catahoula 

Morehouse 

Natchitoches 

St.  John  the  Baptist 

St.  Tammamy 

West  Baton  Rouge 

Winn 


Maryland 
Queen  Annes 


Mississippi 

Chickasaw 

George 

Jackson 

Lamar 

Leflore 

Marshall 


Pearl  River 
Simpson 
Tunica 
Union 
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Missouri 


North  Carolina 


Pemiscott 


Beaufort 

Carteret 

Dare 

Harnett 

New  Hanover 

Northampton 

Perquimans 

Wilson 


South  Carolina 


Tennessee 


Clarendon 
Marion 
Berkeley 
Hanpton 


Carroll 
Henry 


Texas 


Virqinia 


Camp 

Chambers  (alt.  for  Fort  Bend  TX) 
Hardin 

Jasper  (alt.  for  Montgomery  TX) 
Liberty 

Nacogodoches  (alt.  for  Polk  TX) 
Rusk  (alt.  for  San  Jacinto  TX) 
San  Patricio 
Tyler 


Nansemond  ( Suf  f o lk ) 
Northhampton  (alt.  for 

Chesapeake  VA) 
Prince  George 
Surry 
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Central  Great  Plains  Winter 
Wheat  &  Range  Region  (45  counties) 


Kansas 

Chautauqua 

Cloud 

Ford 

Graham 

Hamilton 

Harper 

Kingment  (alt.  for  Canadian  OK) 

Kiowa 

Logan 

Lyon 

McPherson 

Phillips 

Saline 

Sheridan 

Trego 

Wabaunsee 

Wichita 


Nebraska 

Butler 

Hall 

Harlan 

Kearney 

Keith 

Seward 

New  Mexico 

Curry 

Oklahoma 

Beaver 
Caddo 

Cotton  (alt.  for  Ellis  OK) 
Harmon  (alt.  for  Kiowa  OK) 
Jefferson  (alt.  for  McClain  OK) 
Tillman  (alt.  for  Woodward  OK) 


Texas 

Callahan 

Carson 

Cottle 

Dawson 

Garza 

Hale 

Hemphill 

Hutchinson 

Jack 

Lipscomb 
Midland 
Potter 
Scurry 

Stonewall  (alt.  for  Tom  Green  TX) 
Yoakum 
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Data  Center  in  Sioux  Falls,  South  Dakota.    The  dates  of  the  images  were 
determined  by  the  season  of  the  year  in  which  the  changes  would  be  most 
obvious.    In  the  Southeast,  spring  imagery  was  selected  to  show  the 
best  contrast  between  forest  and  cleared  land.    It  was  assumed  that 
in  the  spring,  cropland  would  have  a  bare-aground  response  on  the 
imagery,  due  to  recent  plowing  or  planting,  while  pasture  would  have 
a  slight  vegetative  response.    In  practice  it  was  impossible  to  dis- 
tinguish cropland  from  pasture,  so  this  distinction  had  to  be  aban- 
doned. 

In  the  Great  Plains,  late  summer  and  early  autumn  dates  were 
used.  During  the  dry  season  it  was  thought  irrigated  crops  would 
be  taller  and  greener  than  nonirrigated  crops  or  range.  Distinc- 
tions between  range  and  nonirrigated  cropland  would  be  determined 
by  field  patterns  and  vegetative  differences. 

Spring  imagery  for  the  Southeast  was  only  available  for  1973. 
Because  several  scenes  in  the  Southeast  did  not  have  acceptable 
images,  images  from  the  fall  of  1972  or  1973  were  occasionally 
substituted.    This  occured  despite  the  frequency  of  Landsat  coverage. 

Images  were  ordered  from  EROS  in  a  55.8mm  black  and  white  film 
positive  form  (called  a  chip).    Generally,  for  each  scene,  bands  5 
and  7  were  ordered  for  the  early  year  and  band  5  was  ordered  for  the 
later  year.    This  combination  was  determined  through  experimentation 
at  the  beginning  of  the  project  on  three  counties— Beaver  County, 
Oklahoma;  Ford  County,  Kansas;  and  Catahoula  Parish,  Louisiana.  The 
combination  of  bands  5  and  7,  once  enlarged  and  color  enhanced,  pro- 
vided a  greater  color  volume  as  a  base  to  contrast  against  the  image 
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of  the  latter  year.l/    In  the  Southeast,  changes,  could  be  identified 
from  band  5  alone. 

The  images  were  then  enlarged  in  two  stages  to  an  8X10  film 
positive.    This  was  a  scale  change  of  1:3,300,000  on  the  chip  to 
1:370,000  on  the  film  positive.    Highway  maps  of  each  county  were 
ordered  and  enlarged  to  a  standard  scale  of  1:63,360  (1  inch  per 
mile)  to  use  as  a  base  map  for  each  county. 

Diazo  transparencies  of  the  8X10  film  positive  were  prepared 
using  diazo  color  film  in  magenta,  cyan,  and  yellow.  Interpretation 
keys  were  prepared  based  on  work  done  on  the  three  pilot  counties: 
Ford,  Beaver,  and  Catahoula  Parish  (See  Appendix  1).    Areas  of  change 
between  the  two  years  were  identified,  coded,  and  mapped  on  the  base 
map  via  a  projector.    Land  use  change  acreages  were  then  tabulated 
using  a  prepared  grid. 

Analysis  of  Results 
The  data  collected  from  the  Landsat  analysis  are  presented  in 
tables  2  and  3.    A  small  percentage  of  the  land  area  could  not  be 
interpreted  due  to  cloud  cover.    It  is  difficult  to  discuss  the 
reliability  of  this  data  as  location  problems  in  field  checking 
prevented  definitive  conclusions.    Consequently,  no  regional  expan- 
sions are  attempted,-.*  However,  using  the  field  check  data  and  avail- 
able secondary  data  some  tentative  conclusions  can  be  drawn. 

Informal  field  checks  were  conducted  throughout  the  project  by 
1/  For  further  description  of  this  process,  see  (Skaley,  1977). 
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Table  2.    Landsat  data  for  Southeast  regional  group,  1972-76,  in  acres 

Forestland  cleared  Not  interpreted  Total  land  area 


Alabama 

Barbour 

2,438 

16,840 

570,240 

Choctaw 

3,682 

583,040 

Dale 

2,115 

— 

357,440 

Henry 

8,037 

354,560 

Mobile 

7,521 

8,160 

793,536 

Tuscaloosa 

5,103 

853,120 

Wilcox 

6,147 

575,168 

Florida 

Baker 

1,341 

— 

374,144 

Calhoun 

3,407 

358,912 

Dixie 

3,564 

— 

443,136 

Jefferson 

4,028 

387,328 

Liberty 

1,510 

— 

536,704 

Nassau 

8,661 

416,000 

Washington 

1,862 

— 

374,592 

Kentucky 

Calloway 

609 

— 

245,824 

Maryland 

Queen  Annes 

11,417 

27,820 

240,192 

Missouri 

Pemiscott 

3,841 

315,456 

South  Carolina 

Clarendon 

15,567 

383,552 

Marion 

1,878 

311,936 

Berkeley 

12,340 

710,144 

Hampton 

2,212 

359,424 

Texas 

r* 

Camp 

285 

123,008 

Chambers 

259 

394,304 

Hardin 

1,812 

574,080 

Jasper 

6,078 

580,480 

Liberty 

1,038 

24,960 

755,200 

Nacogodoches 

2,897 

577,280 
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Forestland  cleared  Not  interpreted  Total  land  cleared 


Texas 

Rusk 

711 



600,768 

San  Patricio 

2,183 

438,592 

xyxer 

A   9  AC\ 
4  ,Z4U 

/  ,  OUU 

doo, uy  0 

North  Carol  ma 

rXrctUXvJI.  L. 

X  ,  .700 

coo  70/1 

JZO  ,  /  U4 

C  ^  IT  +■  C±Y~C1  +■ 

v-ctl.  Lcl  t:  U 

1  OfiS 

X  ,  U  O.J 

?9  4S^ 

4?9  91 9 

UaL  tr 

X  ,  OOJ 

Harnett 

1,689 

— 

385,664 

New  Hanover 

230 

118,656 

North  Hanpton 

1,672 

343,040 

Perquimans 

902 

— 

157,696 

wixson 

/uu 

9  9Q  7/1/1 

z jy , / 44 

Tennessee 

Carroll 

265 

— 

381,312 

nci 

74? 

jDj , UU O 

V XLh  ill J.d 

Nansemond  (Suffolk) 

1,989 

* 

Northampton 

446 



140,864 

Prince  George 

1,500 

— 

182,400 

Surry 

1  SQQ 

177    A O  Q 

x  /  /  , 4UO 

LAroI  let 

x  ,  jzu 

47f)  Ql 9 

Cleveland 

541 

27,840 

384,640 

Lee 

1,697 

388,992 

Little  River 

150 

311,040 

Phillips 

1,839 

— 

450,565 

Dy-o  -j  y  -j 

4  41  1 
*± ,  4xx 

499  848 

4ZZ , 040 

Otrv->i.y  Xd 

ui  cu  l  u  xtsy 

4D9 

U  ,  4UZ 

Z  OD , U  OU 

Brooks 

7,796 

313,920 

Clay 

3,598 

128,000 

Clinch 

7,577 

509,760 

Dodge 

3,588 

318,976 

Evans 

1,051 

118,848 

*  No  data  available 
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Forestland  cleared  Not  interpreted  Total  land  area 


Georgia 

Houston 

2,660 

— 

243,200 

Jeff  Davis 

2,956 

— 

211,520 

Macon 

4,474 

— 

257,664 

Pi  £*y~c^ 

C  J.  —  i-  V^v— 

2  157 

218  880 

Pulaski 

2,359 

— 

162,112 

Q^ttii  nol  (=> 

458 

157  184 

Terrell 

3,789 

— 

210,560 

Turner 

2,766 

— 

187,328 

Louisiana 

Beauregard 

3,988 

— 

757,568 

Calcasieu 

3,071 

6,360 

707,200 

Caldwell 

1,451 

— 

352,448 

Cameron 

513 

99,230 

922,368 

Catahoula 

4,973 

— 

474,944 

Morehouse 

9,949 

— 

514,560 

11  554 

828  800 

St.  John  the  Baptist 

201 

— 

145,024 

St«  Tammany 

521 

567,360 

W.  Baton  Rouge 

539 

— 

129,856 

Winn 

9,561 

— 

608,256 

Mississippi 

Chickasaw 

1,689 

— 

323,776 

George 

4,596 

1,118 

307,520 

Jackson 

5,685 

2,534 

470,848 

Lamar 

1,082 

— 

320,000 

Leflore 

2,379 

9,900 

379,008 

Marshall 

4,722 

454,208 

Pearl  River 

624 

529,920 

Simpson 

1,648 

8,640 

375,686 

Tunica 

2,646 

293,056 

Union 

434 

269,952 
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Table  3.    Landsat  data  for  Great  Plains  regional  group,  1972-76  in  areas. 


Range land  to 

Range land  to 

Nonirrigated 

Not 

Total  land  a 

Nonirrigated 

Irrigated 

Cropland  to 

Interpreted 

Jy~y     Pvonl  ^nH 

■LI.  L  •      vl  v_/^_L  CU  H_i 

Kansas 

Chautauqua 

Q89 

41  A  OP.P 
4-L'±  ,  UOU 

Cloud 

Q80 
_7  ou 

AC.C.  r)Af) 
'iDD  /  U4U 

Ford 

1  768 

4,418 

17,220 

IS  824 

698  176 

Graham 

1  007 
J.  ,  UU  / 

1  764 

^  Q  PI 

C70    1 7 A 
D  /  U  ,  ±  /  O 

Hamilton 

i  o  S44 

J-U  ,  D'iH 

6  189 

9  81  4 

614  880 

Harper 

1  09R 

Q94 

1  4  400 
-L*4  ,  4UU 

DJ.Z  ,  D  4u 

Kingman 

JOJ 

1  Q94 

X  ,  Z7  Zfi 

ceo  qcn 

Kiowa 

9  7S9 

d  81  n 

t  ,  OJ.U 

1  1  7S1 

46D  8DD 

Logan 

1    1  96 

40 

*±U 

4  Q7Q 

686  790 
U  oo , / zu 

Lyon 

Z4Z 

D  ,  UJ-D 

COO  9/in 
jjo , Z4U 

McPherson 

77 

S   1  14 

4  140 
*±  ,  J  **u 

S71  11 9 

Phillips 

/  4:? 

9  974 

£  0  9Q 

R74  90S 

J  /  'i/ZUO 

Saline 

355 

332 

313 

— 

460,800 

Sheridan 

460 

3,379 

20,082 

— 

571,456 

Trego 

1,367 

588 

1,767 

576,640 

Wabaunsee 

732 

610 

3,384 

— 

507,136 

Wichita 

1    7  Ql 

91  70£ 
Z± ,  / UD 

490  484 

1ZU , 1u4 

Nebraska 

Butler 

-379  AQf) 

Hall 



973 

8,919 

— 

343,360 

Harlan 

— 

9,655 

355,904 

Kearney 

7,222 

10,254 

327,680 

Keith 

81 

5,355 

8,617 

660,672 

Seward 

9,662 

13,440 

365,632 

New  Mexico 

A     C  /1Q 

4,b4y 

z,4Jz 

op  QA£ 
Zo , 040 

897  792 

uiuanoma 

leaver 

4,075 

484 

o,  ojy 

i  9  nnn 

J.Z  ,  UUU 

1   1  SI  710 

JL,-LJO,  /  OU 

Caddo 

Q4Q 

1 .624 

13,900 

814,080 

Cotton 

439 

V 

305 

416,320 

Harmon 

483 

773 

34,030 

348,800 

Jefferson 

5,571 

328 

2,400 

499,200 

Tillman 

530 

379 

2,894 

576,832 
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Rangeland  to     Rangeland  to 
Nonirrigated  Irrigated 
Cropland  Cropland 


Texas 

Callahan 

901 

75 

Carson 

— 

213 

Cottle 

1,068 

— 

Dawson 

472 

44 

Garza 

529 

— 

Hale 

— 

1,381 

Hemphill 

601 

173 

Hutchinson 

527 

1,207 

Jack 

— 

— 

Lipscomb 

5,180 

6,199 

Midland 

751 

52 

Potter 

Scurry 

374 

Stonewall 

839 

Yoakum 

6,969 

1,903 

Nonirrigated  Not  Total  land  area 

Cropland  to  Interpreted 
Irr.  Cropland 

300  —  548,047 

24,891  —  576,000 

—  —  576,064 
5,215  —  577,280 

88  —  584,896 

22,696  —  626,560 

1,513  —  578,304 

16,100  —  559,936 

—  —  604,928 
4,603  —  597,760 

749  —  600,768 

1,350  —  574,720 

46  —  578,304 

283  —  592,768 

9,763  —  531,200 
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RIL  and  USDA  personnel.    Gounty  Soil  Conservation  Service  (SCS)  per- 
sonnel were  also  consulted  on  an  informal  basis  during  the  interpre- 
tation phase  of  the  project.    However,  a  formal  field  check,  using 
SCS  personnel,  was  attempted  after  the  interpretation  phase  of  the 
project  had  been  completed.    A  10  percent  sample  of  land  use  change 
areas  in  a  sample  of  counties  was  selected  for  ground  checks.  Iden- 
tification maps  were  prepared  with  numbered  sample  sites  for  ground 
checks  by  SCS  county  personnel.    The  SCS  Conservationists  were  asked 
to  determine  the  land  use  for  each  sample  site  in  1972  and  again  in 
in  1976.    A  sample  of  areas  with  no  change  were  included  as  a  check 
on  omissions  in  interpretation.    An  example  of  a  field  check  form  is 
included  in  Appendix  2. 

For  the  6  pilot  counties  field  checked  in  this  manner,  the  re- 
sults were  disappointing. 2/    In  Dare  County,  North  Carolina  where 
logging  was  underway,  clear  cutting  as  a  part  of  a  timber  operation 
was  indistinguishable  from  land  cleared  for  crops.    The  results  from 
the  remaining  five  counties  indicated  that  SCS  field  observations 
and  the  satellite  interpretations  agreed  in  39  of  131  sample  sites, 
for  a  30  percent  accuracy  rate.    Dividing  the  sample  sites  into  change 
and  no  change  sites,  satellite  interpretation  and  field  checks  agreed 
in  21  out  of  32  no  change  sites  (65%)  and  18  of  99  change  sites  (18%). 

Two  explanations  are  possible  for  the  apparent  differences  in 

accuracy  between  change  and  no  change  sample  sites.    One  is  that  map- 

ping  location  problems  prevented  field  checkers  from  locating  the 

2/  The  six  counties  field  checked  in  this  manner  were  Dare,  North 
Carolina;  Cateret,  North  Carolina;  Queen  Annes,  Maryland;  Sheridan, 
Kansas;  Clinch,  Georgia;  and  Union,  Mississippi. 
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exact  site.    Since  the  no  change  sites  were  generally  larger  than  t 
sample  change  sites,  this  would  account  for  the  more  accurate  ident 
ficaticn  of  no  change  sites.    Another  possible  explanation  is  that 
much  of  the  land  use  change  detected  via  satellite  imagery  was  a 
signature  problem  involving  cloud  conditions,  water,  or  crop  con- 
ditions. 

The  possibility  that  Landsat  Imagery  was  not  aligned  with 
county  maps  stems  in  part  from  distortion  in  the  imagery.  Geo- 
graphic correction  is  necessary  on  satellite  imagery  as  the  only 
undistorted  part  of  each  scene  is  the  center,  the  closer  to  the 
edge  of  the  scene  the  greater  the  angular  distortion.    In  addi- 
tion, the  land  use  change  sites  in  this  project  were  particularly 
difficult  to  relate  to  a  map  as  they  were  scattered  and  irregularly 
shaped. 

However,  the  formal  and  informal  checks  indicated  that  there 
were  many  deviations  from  the  signatures  established  by  the  inter- 
pretation keys.    Vegetative  response  not  only  differed  within  the 
region,  but  within  counties,  between  scenes  and  within  the  same 
scene.    Sane  of  the  interpretation  problems  encountered  are  de- 
scribed below: 

a.  Inability  to  distinguish  between  pasture  and  cropland  in 
the  Southeast. 

b.  Fields  of  soybeans  may  appear  as  forested  areas  on  images 
made  late  in  the  spring.    Hie  crop  at  that  time  may  be 
dense  enough  to  be  misinterpreted  as  brush  or  forest. 

c.  Flood  conditions  in  several  counties  in  the  Southeast 
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caused  misinterpretations  of  land  uses. 

d.  Deciduous  forests  in  the  Southeast  may  show  up  on 
October  imagery  as  cleared  areas,  if  the  images  were 
made  after  the  autumn  leaf  fall. 

e.  In  the  Great  Plains,  natural  vegetation  may  have  been 
misclassif ied  as  newly  irrigated  land.    Wetter  soil 
conditions  in  the  latter  year  due  to  recent  rainstorms 
or  a  wet  year  could  cause  natural  vegetation  to  have  a 
similar  spectral  response  as  irrigated  crops. 

f.  Conversions  of  range  to  nonirrigated  cropland  were 
extremely  difficult  to  distinguish  from  the  imagery. 

In  some  instances,  field  boundary  changes  were  the  only 
means  of  identifying  the  change. 

g.  The  use  of  different  fall  dates  may  have  caused  some  mis- 
interpretations of  irrigated  cropland.    If  the  image  of 
the  early  year  is  taken  after  the  harvest  of  the  irrigated 
crop,  the  land  may  appear  as  nonirrigated.    If  the  image 
of  the  later  year  is  taken  earlier  in  the  season,  before 
harvest,  the  same  area  would  appear  as  irrigated  and  inter- 
preted as  a  land  use  change. 

Secondary  data  might  help  evaluate  the  project  results.  However, 
no  secondary  county  data  is  available  on  clearing  in  the  Southeast  for 
the  1972-1976  time  period.    Some  secondary  data  is  available  on  irri- 
gation in  the  Great  Plains  states  which  corresponds  generally  to  the 
1972-76  time  period  (table  4).    This  data  was  collected  from  county 
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Table.  4    Irrigation  changes  in  selected  counties  recorded 
by  irrigation  survey  and  Landsat  data. 


Change  in  Irrigation  Change  in  Irrigation 

Irrigation  Survey  Landsat 
 1972-75  1972-76  

Kansas 

Chatauqua  -500  — 

Cloud  1,000  980 

Ford  26,500  21,638 

Graham  9,300  7,745 

Hamilton  4,520  16,196 

Harper  415  924 

Kingman  4,926  2,509 

Kiowa  2,356  18,563 

Logan  2,700  5,019 

Lyon  -160  5,015 

McPherson  14,720  5,134 

Phillips  788  8,403 

Saline  2,896  645 

Sheridan  42,968  23,461 

Trego  470  2,355 

Wabaunsee  -170  3,994 

Wichita  -35,040  23,487 

Total  77,689  146,068 


Texas 

Carson  5,875  25,104 

Dawson  -38,500  5,259 

Garza  -815  88 

Hale  34,220  24,077 

Hemphil  6,500  1,686 

Hutchinson  22,000  -  17,307 

Lipscomb  13,200  10,802 

Potter  1,700  •  1,350 

Scurry  1,000  46 

Yoakum  v28,000  11,666 

Total  73,180  97,385 
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Irrigation  Survey 
1972-76 


Landsat 
1972-76 


Oklahoma 

Beaver 

9,000 

6, 823 

Caddo 

52,500 

2,573 

Cotton 

-1,697 

305 

Harmon 

560 

773 

Jefferson 

100 

328 

Tillman 

8,800 

3,273 

Total 

69,263 

14,075 
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extension  agents  for  the  most  part  and  represents  informed  esti- 
mates of  irrigation  activity,  not  field  surveys.    In  addition,  this 
data  represents  net  irrigation  change,  while  the  Landsat  data  encom- 
passes only  additions  to  irrigated  land.    In  some  areas,  substantial 
acreages  were  removed  from  irrigation  during  the  1969-74  time  period. 
Even  in  areas  with  no  recorded  irrigation  declines,  conversions  to 
irrigation  may  have  been  offset  by  withdrawal  of  irrigation  systems 
on  other  cropland. 

There  is  considerable  variation  between  the  Landsat  data  and 
the  irrigation  survey  data.    These  variations  are  difficult  to  eval- 
uate, however,  due  to  the  different  techniques  used  and  the  differ- 
ence between  net  and  gross  data.    Interpretation  of  both  field  check 
and  secondary  data  would  appear  to  indicate  that  these  are  interpre- 
tation problems  in  the  Landsat  data.    Even  if  the  hypothesis  that 
map  location  problems  were  the  reason  land  use  change  areas  could 
not  be  verified  in  the  field,  then  measuring  land  use  change  via  Land- 
sat only  achieved  65  percent  accuracy. 

Guidelines  for  Future  Research 
There  are  certain  limiting  factors  discovered  in  the  course  of 
this  study  that  should  be  considered  for  any  future  work.    Some  of 
these  limitations  are  inherent  in  the  Landsat  technology,  while 
others  may  be  improved  upon  based  on  the  experience  of  this  study. 

The  first  limitation  encountered  in  this  study  was  the  problem 
of  cloud  cover.  This  was  a  problem  in  both  of  the  regions  studied, 
although  it  was  more  serious  in  the  Southeast.    Because  of  cloud 
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cover  alternate  counties  had  to  be  selected,  departing  from  the 
original  sample  and  introducing  some  bias  if  the  data  were  to  be 
expanded  to  a  regional  basis.  This  problem  was  compounded  through 
the  multi-period  approach  of  the  project.    Comparable  dates  were 
often  not  cloud  free  for  both  1972  and  1976.    This  required  substi- 
tutions, to  the  extent  of  spring  for  fall  imagery  in  some  cases, 
and  prevented  complete  comparability  in  the  time  periods  used. 
Interpretation  problems  were  thus  compounded.    Even  with  the  use  of 
alternate  counties,  a  small  percentage  of  the  land  area  could  not 
be  interpreted  because  of  cloud  cover. 

Cloud  conditions  are  one  factor  inherent  to  any  use  of  Landsat 
imagery.    Any  multi-period  approach  to  a  problem  using  Landsat  will 
be  seriously  limited  if  cloud  cover  prevents  comparable  time  periods 
from  being  analyzed.    Furthermore,  any  sampling  will  be  difficult 
using  Landsat  judging  from  the  results  of  this  study  and  the  use  of 
alternate  counties. 

Another  limiting  factor  in  this  study  was  the  use  of  one  inter- 
pretation key  for  all  the  counties  in  a  region.    In  actuality,  field 
checking  showed  that  there  were  many  deviations  from  the  signatures 
established  by  this  key. 

These  variations  in  interpretation  can  be  reduced  by  careful 
field  checking  at  various  stages  in  the  project.    The  selection  of 
imagery  should  be  checked  first  to  make  sure  no  unusual  weather  con- 
ditions such  as  flood  or  drought  were  present  that  might  result  in 
misinterpretations. 

Second,  interpretation  keys  should  be  based  on  a  knowledge  of 
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known  sample  sites.    A  variety  of  these  sites  should  be  located 
throughout  the  study  region.    Mid-project  field  checks  in  other  se- 
lected sites  should  be  conducted  to  correct  for  any  unforseen  inter- 
pretation problems.    The  telephone  interviews  conducted  in  this 
study  were  not  sufficient.     It  is  doubtful,  however,  that  all  inter- 
pretation problems  can  be  eliminated  so  a  final  field  check  should 
be  conducted  to  assess  the  accuracy  of  the  interpretations. 

The  failure  of  the  field  checks  in  this  study  to  adequately 
assess  the  project  stemmed  largely  from  the  'uncertain  alignment  of 
Landsat  imagery  with  county  maps.  Procedures  for  registration  of 
scenes  to  map  base  have  been  developed  by  the  Statistics  section 
of  the  USEA  for  use  in  their  digital  programs.    New  Landsat  techno- 
logy will  provide  data  that  is  already  corrected  for  geographic  dis- 
tortion, making  registration  to  a  map  base  simpler. 

In  addition,  the  land  use  change  sites  in  this  project  were  dif- 
ficult to  field  check  because  they  were  scattered,  irregularly  shaped, 
and  not  always  visible  from  a  road.    The  Statistics  project  uses  as 
field  sites  known  sample  plots  from  their  June  Enumerative  Survey  to 
establish  a  ground-truth  data.  base.    In  this  manner,  Landsat  imagery 
is  tied  to  key  plots  where  ground  data  is  clearly  defined.    An  ap- 
proach of  this  type  would  appear  to  be  preferable  to  establishing 
ground  data  after  the  interpretation  process  is  completed. 

Field  checking  procedures  are  an  important  component  of  any  Land- 
sat project.    The  difficulties  encountered  in  this  study  can  probably 
be  reduced  in  future  projects. 
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Conclusion 


The  results  of  this  study  are  somewhat  ambiguous  as  the  field 
checking  process  did  not  allow  definitive  conclusions  as  to  accuracy 
of  the  data.    However,  it  appears  that  there  are  significant  problems 
in  interpreting  satellite  imagery,  among  them,  unusual  weather  condi- 
tions, difficulty  in  generalizing  interpretations  over  large  areas, 
and  confusion  between  spectral  responses.    In  addition,  the  problem 
of  cloud  cover  makes  any  sampling  extremely  difficult.    In  this  study, 
13  percent  of  the  counties  in  the  Great  Plains  and  19  percent  in  the 
Southeast  were  alternates,  introducing  some  sample  bias.    In  this  study, 
the  use  of  two  time  periods  compounded  the  cloud  cover  and  interpreta- 
tion problems.    Nevertheless,  Landsat  can  offer  unique  information  for 
the  agricultural  sector.    Combining  Landsat  data  with  other  measurement 
techniques  may  be  the  most  promising  use  of  satellite  imagery,  but 
where  no  other  information  is  available  satellite  imagery  alone  could 
be  useful,  if  qualified  by  accuracy.    In  both  cases,  careful  field 
checks  should  be  an  integral  part  of  the  interpretation  process.  It 
is  hoped  that  the  experiences  of  this  study  will  enable  researchers 
in  future  projects  to  make  maximum  use  of  the  information  potential  of 
satellite  imagery. 
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Appendix  1.    Identification  keys  for  Great 

Plains  and  Southeastern  Counties 
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INTERPRETATION  KEY  FOR  WESTERN  COUNTIES 


Land  Use;        Irrigated  Cropland 
Interp.  Pivot  Systems  Other 


Non- Irrigated 
Cropland 


Rang el and 


Physical  ' 
Characteristics 


Dark,  solid 
color  distinct 
circles;  size: 
1/2"  in  diara. 
at  1:63,360 


Dark,  solid 
•color  sharply 
defined  boun- 
daries--square 
or  rectangular 
size:  1/2" 
square  or 
smaller  (unless 
neighboring 
sections  are 
irrigated) . 


Rectilinear 
patterns 
regular 
boundaries 


General ly 
irregular 
boundaries,  but 
may  occur  in 
patches  in  cropped 
areas,  generally 
large  areas, 
mottled  appearance 
drainage 'channels ' 
may  be  visible, 
but  faint. 


Color 

Signatures 
Cyan- 5 
Yellow-7 


Late  Year 
Magenta- 5 


Composite 


Dark  blue 


Dark  blue 


Dark  magenta 


new  irrigation 
dark  magenta 
old  inrigation- 
dark  purple 
discontinued 
irrigation- 
dark  cyan 


Dark  magenta 


Tone  range 
from  pale 
yellow  to 
yellow-green 
to  light 
green 

tone  range 
from  pale 
pink  to  med 
pink 

old:  light 
yellow  brown 
to  light 
green  brown 
new:  see  note 
below 


Dark  green- 
maybe  slightly 
mottled. 


medium  dark 
magenta,  maybe 
slightly  mottled 


old:  dark  green- 
brown,  maybe 
slightly  mottled 


NOTES : 

Conversions:    Locate  new  irrigation/cropland  given  above  guidelines  and 
identify  prior  use  based  on  diazos  of  early  year. 

Rangeland  to  Non-Irrigated  Cropland:    Shows  a  very  subtle  visual  change — 

distinctions  may  appear  primarily  as  a  boundary  shift  evident  on 
late  year  diazo  as  xompared  to  early  year.    Would  occur  primarily 
along  the  interface  between  rangeland  and  cropland  (evident  in  diazo 
of  early  year).    The  area  will  be  lighter  on  the  later  year  diazo, 
thus  the  dark  green  of  the  early  year  will  be  visible. 
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NOTES:  (continued) 


Irrigation  equipment  can  be  moved — so  may  see  a  system  which  has  gone 
out  of  use  (cyan)  near  a  new  one  (magenta) .    This  happens 
occasionally.    We  will  make  no  adjustments  in  the  tabulation 
for  this. 

Imagery:    In  some  cases,  Band  5  was  ordered  for  early  year  and  Bands  5 
arid  7  for  late  year.     Interpretation  is  then  reversed--new 
irrigation  is  cyan,  not  magenta  and  so  on. 
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INTERPRETATION  KEY  FOR  SOUTHEASTERN  COUNTIES 


Land  Use; 


tlterp 


Forest 


Cropland 


Physical 
Characteristics 


Color  Signatures 
Early  Year 
Cyan-Band  5 


Late  Year 
Magenta-Band  5 


ComDosite 


Dark  areas,  no  regular 
boundaries,  large, 
may  follow  stream/river 
patterns,  color  may 
show  some  mottling 


Dark  blue,  fairly  solid 
color,  but  may  show 
some  mottling 

Dark  magenta,  fairly 
solid  color,  but  may 
show  seme  mottling 

Forest  both  years — 
dark  purple 


Rectilinear  areas, 
regular  boundaries 
may  be  small  areas 


lear  to  light  blue 


Clear  to  light  pink 


Cleared  both  years — 
clear  to  light  purple 


NOTES : 


Conversions:  Areas  forested  in  early  year,  cleared  in  late  year 

visible  as  dark  blue  (early  year  forested  response)  areas 
showing  through  a  clear- or  pale  magenta  area  (non-forested 
response ) . 

Colors:     In  some  instances,  the  early  year  Band  5  is  in  magenta 
(instead  of  cyan)  and  the  late  year  Band  5  is  in  cyan 
(not  magenta).     In  such  cases,  they  appear  as  magenta 
areas . 

Band  7  was  ordered  for  all  areas,  but  the  yellow  diazo  of  this 
band  showing  moisture  responses  often  conceals  detail 
and  makes  the  Interpret at ion  more  difficult.     For  this 
reason,  3and  7  (in  yellow)  was  used  only  where  it  does 
not  interfere  with  interpretation. 

Flooding:     Several  areas  in  the  study  may  have  flood  conditions 
evident  on  the  image  (obvious  from  river  boundaries  and 
widths).     This  may  reduce  the  accuracy  of  the  interpret  at  io] 
since  wet  conditions  can  conceal  detail. 

Clearcutting  is  a  common  forestry  management  practice  in  pine 
stands  in  this  region.     Clearcut  areas  which  have  not 
been  replanted  would  be  interpreted  as  areas  cleared 
for  cultivation.     SCS  contact  people  would  be  able  to 
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verify  if  clear-cutting  is  done,  but  may  not  be  able 
to  identify  locations. 

Land  which  has  been  cleared  for  housing  developments  would 

likely  be  interpreted  as  areas  cleared  for  cultivation. 
As  with  clear-cut  forest  lands,  we  would  nor  be  able  to 
'   identify  the  intended  use  of  a  given  cleared  area. 
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Appendix  2.    Sample  Field  Verification  Form 
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FIELD  VERIFICATION  TALLY  SHEET 
VSDA/ZSCS  IANDSAT  SATELLITE  AGRICULTURAL  I 
CHANGE  STUDY 


LAND  USE 


County  CVHt  *.t 


State  ^/  / 


District  Conservationist  fl  .  H,,,,.., .  ,l     £   ^  ^ 

Office  Address  2^^^  'P,..\A^      tr       .  ~ 


QUESTIONS : 

1    If  the  site  was  i 


was  i„  crops,  what  type  of  crop  was  grown?        A  M  „  U    C-  ,  A 


2.    If  the  site  was  classified  as  forest,  was  cl. 
the  time  periods?    Please  specify. 


/ 


arcutting  involved  in  either  of 


3"  *2TTr  SPeCial  ClrCUmStanC6S  (fl0Ods'  that  preve„ted  . 

norn*!  land  use  pattern  durin,  either  of  the  tin*  Periods,  M.J, 


e  specify. 


Additional  coinnents  about  site: 


i 
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